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CONCLUSIONS 


A difference equation has been derived that describes the performance of a 
digital phase-locked loop of order one or two. A computer program has been written 
that simulates the difference equation and hence simulates the DPLL itself. Using this 
simulation program, the performance characteristics of the DPLL for an input signal 
corrupted by noise with respect to several loop parameters were determined. In 
particular, the steady-state phase error as a function of signal-to-noise ratio for a 
constant phase input and the mean time-to-lock for a initial phase offset have been 
determined by means of this simulation program. 

II. INTRODUCTION 

The lack of suitable mathematical models for digital phase-locked loops of 
order higher than one makes it necessary to evaluate such loops solely by means of 
simulation. Two notable derivations by Reddy and Gupta [1] and Gill and Gupta [21 
of simulation models have appeared in the literature recently that describe the simulation 
of first-and second-order digital phase-locked loops of a somewhat general configuration. 

The purpose of this paper is to apply the results of [1] and [2] to a particular DPLL configuration. 
Slight changes to ihe material referenced will be required to accommodate the particular DPLL 
configuration considered here; therefore, the simulation equations used in this paper will 
be derived in detail. 

In Section IV these simulation equations will be applied to various loops to 
determine the effects of the several loop parameters upon the output phase. Specifically, 
the steady-state phase error and the mean time-to-lock for various loops will be examined. 

III. SIMULATION EQUATIONS 

The DPLL used in the simulation is given in Figure 1 . For this loop, the input 
is sampled by the phase detector at the positive-going zero crossing of the reference 
clock. The output of the phase detector is quantized to plus or minus one depending on 
the sign of the sample. Sampling in this manner, the phase difference between the reference 
clock and the input signal is essentially reduced to a determination of phase lead or lag 
of the reference clock with respect to the input. The two forward paths represent a filter 
for the output of the phase detector in the sense that they alter the detector outputs to 
create a signal suitable for adjustment of the reference clock phase. 

Operation of the loop can best be described by considering the special case in 
which the noise is assumed to be zero. The input signal then is 

s(t) = sinu t 
o 

let the reference clock be quantized to n distinct phase states so that the reference clock 
can be described by 
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r(t) = sin (« Q t - ~ +tt ) i-l,2,...N 

assume that initially the reference clock is given by 

r(t) - sin (u Q t - tt ) 

the reference clock then lags the input by 


it Radians. 

2N 

The relative phase of the two signals is depicted in Figure 2. Since the input is sampled 
at the positive going zero crossing of the reference clock, from Figure 2, the first sample 
is positive in sign and the output of the phase detector is +1 . For the forward path of 
Figure 1 containing the A 1 block, this +1 output is multiplied by the value A 1 radians. 
The phase detector outputs are also summed in the other forward path, multiplied by the 
constant A 2, and then passed through the nonlinearity. For the case considered here, 
after the first sample is taken the value at the output of the nonlinearity is A 2 radians. 
The signal Y(k) is then used to alter the phase of the reference clock by 

A 1 + A 2 Red. 

before the next sample is taken. Obviously, A 1 and A 2 must be related to the 
quantization of the loop. If at the second sample, the reference clock still lags the 
input the reference clock phase will be increased by 

A 1 + 2 *A2 Rad. 

This continues until the { th sample, at which time the reference clock phase leads the 
input signal phase and the phase is altered by 


- A1 + A 2 I x(i) 
i=l 


Finally as the loop achieves lock, the value of the output of the nonlinearity is less 
than or equal to A 2 for all successive samples. 

Following the derivation of (2), the k^ sample occurs at time t(k) and the time 
between samples is 

T(k) = t(k)-t(k-1) (1 

Since the output of the phase detector occurs only at the times t(k), then X(k) is 
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( 2 ) 


X(k)- sgn {A sin C« t+9(t(k))] + n [t(k)l } 
c o 

Also, let the k™ sample interval be given by 

T(k) = T - Y(k) (3) 

where T is the period of the input signal. Then the actual time of the k^ sample is 


t(k)= l T(k) 

i-l 

k 

= kT - I y (i) 
i=0 


and the k^ phase detector output is 


k-1 


X(k) - sgn {A sin [9(t(k))-u I Y(i)] + n Ct(k)] } 
c Oj-Q 


If the phase error is defined as 


k-1 


0[t(k)]= 9 Ct(k)]-<^ I Y(i) 

i- 


° 5=0 


then 


where 


<t>(k+l) - 4>(k) = 9(k+l) - 9(k) -u Y(k) 


9 Et(k)3 = 9(k) and 0 Et(k)3 = 0(k) 

But Y(k) is defined as 

k 

Y(k) = A ,X(k) + f { A I X(i) } 

1 1 i=0 

so that Equation (7) can be written as 

0(k+l) - 0(k) = 9(k+l) - 9(k) - A ^ sgn [4>(k) + n(k)] 


(4) 


(5) 


( 6 ) 


( 7 ) 


( 8 ) 


- u f { A 0 I sgn [0(1) + n(i) ] } 
o l ;=0 


(9) 


This final equation then provides a means of simulating the phase error of the loop, and 
hence the phase output of the loop, for any function of phase input. 
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IV. DPLL SIMULATION 


Using Equation (9), the DPLL of Figure 1 was simulated via a computer program 
to determine the characteristics of the DPLL for a signal plus wideband Gaussian noise 
input. In particular, the steady-state phase error was determined for a constant input 
phase of zero radians and the mean time-to-lock for an initial phase offset were deter- 
mined. The programs used for the simulation may be found in the Appendix. 

Figures 3, 4, and 5 give the steady-state phase error response for different loop 
constants. In Figure 3, the reference clock phase was quantized to 32 distinct values 
while the constants Al and A2 were both set to 2tt/128 radians, where 128 is the number 
of distinct loop states. The value of A1 (the nonlinearity) was varied to values of 0, A2 
and 2*A2. For the case of Al = 0, the loop was acting as a first-order loop only, while 
for Al = A2 and 2*A2, the loop is capable of second-order phase updating. From Figure 3, 
it is apparent that for signal-to-noise ratios less than 20 dB that significant steady-state error 
is introduced by the second-order updating of the loop. For signal-to-noise ratios higher 
than 20 dB, the performance of the loop in terms of s*oady-state error is essentially the same 
for all three values of Al . The limiting value of the steady-state phase error for signal-to- 
noise ratios higher than 20 dB is a function of the number of quantized values of the loop 
output phase. As the signal-to-noise ratio decreases below -40 dB, the loop output error 
approaches that as would be found with a uniform phase distribution regardless of the value 
of Al . Figure 4 is a similar plot except that Al and A2 have been set to 2tt/256 radians. 
Comparison of Figures 3 and 4 reveals that for corresponding values of Al and SNR, the 
steady-state phase error given in Figure 4 will be less than or equal to that given in Figure 
3. This is caused by the smaller values used for Al and A2 in the case given in Figure 4 
which essentially indicates that the phase detector outputs are filtered more for this case. 
Figure 5 gives the steady-state phase error for a loop in which the output phase is quantized 
to 64 distinct values and the constants Al and A2 are set to 2tt/512 radians. Notice that 
in this case, for signal-to-noise ratios above 20 dB the steady-state error approaches a value 
one-half that of Figures 3 and 4 and is less than that obtained for Figures 2 and 3 for all 
corresponding values of Al and SNR. 

The mean time-to-lock for the DPLL's considered above were also determined for 
a signal-to-noise ratio of 30 dB to check the settling time of the loops. In Figure 6, the 
mean time-to-lock is given for the same loop constants as were used in Figure 3. Notice 
that for small values of initial phase error that the second-order loop with Al = 2 actually 
takes longer to lock than for the other two cases. This is caused by the maximum value 
building up in the summation block of the second-order forward path and requiring 
additional samples to settle out. This same characteristic is seen in Figure 7 where 
the loop constants are the same as used in Figure 4. In Figure 8, the mean time-to- 
lock for a loop quantized to 64 phase states is given. 
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Figure 3. Standard Deviation of Steady-State Phase Output. 
M = 4, N = 32, = A 2 = 0.049 Rad. 
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• A1 — 0 
A A1 = 1*A 2 
■ Al= 2*A 2 
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Figure 4. Standard Deviation of Steady-State Phase Output. 
M= 8, N = 32, A = A 2 = 0.0245. 
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Mean Time To Lock (Samples) 



Figure 7. DPLL Transient Response. 

M = 8, N = 32, A = A =0.0245 
SNR = 30 dB. ■ 
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APPENDIX 


A. Steady-State Error Program 


DIMENSION TOT (. 64 ) 

D A T A N S E E D v P I y STD y TOT/ 1 37 y 3 ♦ 1 4 159265 rO.» 64*0 * / 
INTEGER STATE 
WRITE* Ay 100). 

100 FORMAT (IX y 'ENTER M-REGISTER LENGTH') 

REAIKSy 101 ) <M) 


10.1. FORMAT (15) 

WR I TE ( A r 1 02 ) 

1 02 FORMAT < 1 X y ' ENTER N-REGISTER LENGTH ' ) 
READ ( 5 y 101 ) <N) 

MN*M#N 
WRITE (6 y A ) 

6 FORMAT (IX y 'ENTER SNR IN DB') 

READ (5 y .1.04 ) (SNR) 

A --SORT (2* )*10.**( SNR/20 ♦ ) 

WRITE (A* 103) 

103 FORMAT ( IX y 'ENTER W*DEL1 ' ) 

READ (5 y 104) (WDEL1 ) 

104 FORMAT (FI 0 ♦ 0 ) 

WRITE* Ay 105) 

1 05 FORMAT ( 1 X y ' ENTER W*DEL2 ' ) 

READ (5 y 104) ( WDEL2 ) 

WRITE (Ay?) 

7 FORMAT ( IX y 'ENTER MAX 2ND ORDER EFFECT') 
READ (5y 104) ( A 1 ) 

PHI“"1 ♦tfPI/MN 
K«PHI/(M*WDEL1 ) 


ORIGINAL PAGE IS 
OE POOR QUALITY? 


PHIQ«FI/N 

if(ph:im..t«o) phio~-p:i:/n 

IF ( K . GT ♦ 0 ) PHIQ*PI/N+K#2 . *P I/N 
I F ( K ♦ I... T ♦ 0 ) P H I Q « - P I / N f 2 . *K * P I / N 
TEMP2*0 . 

IF ( TEMP 2 ♦ EQ ♦ 0 ♦ ) GO TO 5 
TEMP 2 :r - ADS ( TEMP 2 ) / TEMP 2 
5 DO 10 I«.t y 1000000 

CALL GAUSS ( NSEED y 1 . y 0 ♦ y 0 ) 

TEMF > *A#SIN(PHIQ)+V 
IF (TEMP. EG.. 0. ) GO TO 20 
T E M P* A B S ( T E M P ) / T E M P 
2 0 T E M P 2 T E M P 2 1 T E M P 

I F ( TEMP 2 ♦ GT ♦ A 1 ) TEMP2* A 1 
1 F ( TEMP 2 ♦ LT * -A 1 ) TEMP 2™-- A 1 
P H I a P H I w D E L 1 * T E M P - W D E L 2 * T E M P 2 
IF ( PHI ♦ GT ♦ PI ) PHiaPHI-2>*FI 
I F ( P H I . L T ♦ P I ) P H I « P H I + 2 ♦ * P I 
K ™ P H I / ( M * W D E L 1 ) 

PHIO“P.T/N 

I F ( PHI ♦ LT ♦ 0 4 ) PHIQa-PI/N 
I F ( K ♦ GT ♦ 0 ) PI-1 1 QaP I /NT 2 . *K*P I /N 
IF ( K . LT * 0 ) PH I Q-'-PI/N-f 2 ♦ *K*P I/N 
T E M I"' 1 « ( P H 1 G + P I ) / ( 2 . * P I / N ) + 0 ♦ 5 1 
8TD-STD4 ( PHI G ) ##2 
STATE :;:: TEMPI 

TOT ( STATE ) a TOT ( STATE ) + 1 

I F ( MOD ( I y .1 0000 ) * NE . O') GO TO 1 0 

Rial 

STDOaSQRT ( STD/RI ) 

WR I I E (Ay 1 ) ( I y ST DO ) 

1 FORMAT ( 1 X y 7 NO'. OF ITERATIONS^' , IAy5Xy 'PHASE OUTPUT STD DEU~~ ' 
1 »E12. 4) 

10 CONTINUE -14- 


STOP 

END 

SUBROUTINE GAUSS<IX*SjAM*V) 
A~0 ♦ 

DO SO 1-1 ,12 
CALL RANOU(IX»IY»Y> 

IX~IY 
50 A*A+Y 

V~< A-6 t ) XcS+AM 

RETURN 

END 

SUBROUTINE RANDU ( IX 1 1 Y r YFL ) 

IY*IX#65539 

IF < IY ) 5*6,6 

5 IY«IY+2147483647+i 

6 YFL~IY 

YFL--YFL*0 * 46566 1E--9 

RETURN 

END 


R? 



■B. Transient Response Program 


8 


3 F 


7 


20 1 


202 

203 


20 


1ST AND 2ND ORDER DPLL SIMULATION 
DATA N S E E D y P I / 1 37 y 3 ♦ 1 4 1 5 9 2 6 5/ 

DIMENSION TOT < 64) 

INTEGER STATE 
N SEED— 13 9 
WRITE <6v 8) 

FORMAT (1XF ENTER M -REG 1 9 TER LENGTH ' ) 

READ < 5V3 ) < M ) 

WRITE ( 6 y 9 ) 

FORMAT < IX y ' ENTER N-REGISTER LENGTH ' ) 

READ < 5 1 3 ) ( N ) 

ORMAT< 15) 

MN-M*N 
WRITE < 6 y 6 ) 

FORMAT < 1 X y ' ENTER SNR IN DB ' ) 

READ < 5 r.2) < SNR ) 

FORMAT ( FS ♦ 0 ) 

A™ SORT < 2 ♦ ) * 1 0 ♦ ** < SNR/20 * ) 

WDEL1«2.*PI/MN 
WDEL2-2 ♦ TPI/MN 
WRITE (6,7) 

FORMAT ( IX .* ' ENTER MAX 2ND ORDER EFFECT ' ) 

READ < S f 2 ) < Al ) 

WRITE (8*201 > (M>N) 

FORMAT ( ' 1 ' * ' M "-REGISTER LENGTH - ' y 13 y 5X y ' N-REGISTER LENGTH - ' 
1 v 14) 

WR I TE ( 8 y 202 ) ( WDEL 1 y WDEL 2 ) 

FORMAT (1XF W*DEL 1 - ' , El 2 * A r 5X y 7 W*DEL2« '*E12.4> 

WRITERS y20 3) (SNRyAl ) 

FORMAT (1XF SNR- * y F6 . 0 y 5X * ' 2ND ORDER MAX- ' t F5 . 0 > 

PHIM-PI/N 
NN-N/2- I 
DO 200 MM-lvl6 
OAR'-O ♦ 

RMEAN-0 ♦ 

DO 100 J >J— 1 y 200 
PH I -PH I M 
PHI I— PHI 
K*F?HI/<M*WDEL1 ) 

PMIQ-PI/N 

IF ( PHI ♦ I...T ♦ 0 ) PHIQ--PI/N 
IF < K ♦ GT * 0 ) PHIO-PI/NMOK2 ♦ TP I/N 
IF < K ♦ LT * 0 ) PH IQ- -P I /NT 2 * TKTPT/N 
TEMP2-0 ♦ 

I F < TEMP 2 * EG , 0 . ♦ ) GO TO 5 
TEMP2- ADS ( TEMP2 ) /TEMP 2 
DO 10 1-1 y 1000000 
CALL GAUSS ( NSEED y 1 * y 0 * v 0 ) 

TEMP - A T S I N < P H I Q > •}• 0 

I F < TEMP * EG .(),) GO TO 20 

TEMP- ADS < TEMP ) /TEMP 

T E M P 2~ T E MP 2 + T E M P 

IP < TEMP 2 . GT ♦ Al ) TEMP2-A1 

IF < TEMP2 ♦ LT . -Al ) TEMP2»AA1 

PHI -PH I -WDEL 1 TTEMP-WDEL2TTEMP2 

IF< Al ♦ EG . 0*0 ♦ AND ♦ PHI ♦ LE . PI/MN ♦ AND . PHI ♦ GE * -PI/MN > GO TO 11 


****QUU$g 
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K*PHI/(M*WDEL1 ) 

'PHIQ*PI/N 

IF ( PHI ♦ LT * 0 ♦ ) PH I Q-*"P I/N 

IF ( K ♦ 8T * 0 ) PHI Q=P I /N+2 . *K*P I /N 

I F ( K ♦ LT ♦ 0 ) P H I Q * •- P I /N+2 ♦ * K * P I / N 

10 CONTINUE 

11 MAR*VAR+I**2 
WRITE <6 *4) (I) 

4 FORMAT < IX y IS) 

R MEAN- R MEAN + 1 

100 CONTINUE 

RMEAN* < RMEAN/200 ♦ ) 

MAR- ( MAR/200 ♦ ) ~RMEAN**2 
WRITE<8> 13) 

13 FORMAT C IX » //f < MEAN NO. OF SAMPLES TO LOCK <200 TRIALS) / ) 
WRITE <8. 14) (PHIM) 

14 FORMAT < IXr / PHIM= / >F5. 1 * / RAD ' ) 

WRITE ( 8 f 15 ) < RMEAN * MAR ) 

1 5 FORMAT < 1 X , ' SAMPLES * ' » E 1 1 ♦ 3 * ' W I TH A MAR OF ' y E 1 1 * 3 ) 

200 PH'IM*PHIM+2.*PI/N 

STOP 

END 

SUBROUTINE GAUSS ( IX r S y AM r M ) 

A * 0 . 

DO 50 1*1 , 12 
CALL RANDU < IX'y I Y > Y ) 

IX-IY 
50 A*A+Y 

M* < A -6 ♦ )*S+AM 

RETURN 

END 

SUBROUTINE RANDU (IX y I Y r YFL ) 

I Y* IX*65539 
IF(IY) 5 t6,6 

5 I Y « I Y + 21474836 4 7 + 1 

6 YFL-I Y 

YFL*YFL*0 . 465661E-9 

RETURN 

END 
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